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Abstract 

The Canada-France-Hawaii Telescope Legacy Survey, specifically the Very Wide 
segment of data, is used to search for possible main-belt comets. In the first data 
set, 952 separate objects with asteroidal orbits within the main-belt are examined 
using a three-level technique. First, the full-width-half-maximum of each object is 
compared to stars of similar magnitude, to look for evidence of a coma. Second, 
the brightness profiles of each object are compared with three stars of the same 
magnitude, which are nearby on the image to ensure any extended profile is not 
due to imaging variations. Finally, the star profiles are subtracted from the asteroid 
profile and the residuals are compared with the background using an unpaired T- 
test. No objects in this survey show evidence of cometary activity. The second survey 
includes 11438 objects in the main-belt, which are examined visually. One object, 
an unknown comet, is found to show cometary activity. Its motion is consistent with 
being a main-belt asteroid, but the observed arc is too short for a definitive orbit 
calculation. No other body in this survey shows evidence of cometary activity. Upper 
limits of the number of weakly and strongly active main-belt comets are derived to 
be 630 ±77 and 87±28, respectively. These limits are consistent with those expected 
from asteroid collisions. In addition, data extracted from the Canada-France-Hawaii 
Telescope image archive of main-belt comet 176P/LINEAR is presented. 
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1 Introduction 



Historically, asteroids and comets have been thought as two separate pop- 
ulations in the Solar System. Having differing volatile fractions, a common 
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distinction between the two is comets display a coma, while asteroids do not. 
In addition, the two populations differ dynamically. A standard criterion used 
to differentiate between a cometary or asteroidal orbit is the Tisserand pa- 
rameter with respect to Jupiter, 



Tj = ^ + 2 cosJ-(l-e2), (1) 
a V aj 



where aj is the semi-major axis of Jupiter, and a, e, and i are the semi-major 
axis, eccentricity, and inclination of the object, respectively. Comets typically 
have Tj < 3, while asteroids have Tj > 3. 



Recen tly, su rveys such as those d iscussed in Fernandez et al.l (120051 ). IJewitt 



(|2005j), and iLicandro et al.l (120081 ) have discovered a significant number of 
asteroids in comet-like orbits. There are also objects in asteroid- like orbits 
that s how bursts of cometary activity (7968 = 133P/Elst-Pizzaro; lElst et al 



( 19 96J)) o r are a ssociated with meteor streams (near Earth object 3200 Phaethon; 
Whipple! (ll983T )). These observations have made the boundary between comets 
and asteroids less obvious. Intermediate objects may be comets that are ex- 
tinct, dormant, or dead. Conversely, they may be asteroids with higher volatile 
content. 



Objects that display activity but are in asteroid-like orbits are known as main - 
belt comets (MBCs) or activated asteroids (AAs) (IHsieh and Jewittl. 1200 6^. 



Thes e objects are most likely native to the asteroid belt (jFernandez et al 



2002h and may be ac tivated by a collisi o n with a smal l body ( Boehnhardt et al. 



19981 : iTothl . l2000h . IHsieh and Jewittl fl2005l . l2006aT ) carry out the first sur- 



vey to search for MBCs a nd, as of th i s writ ing, only thre e MBCs have been 
found: 133P/Elst-Pi zzaro feist et all Il996h. P/2005 Ul fe.ead et all I2005T ) 
and 176P/LINEAR dHsieh and Jewittl . boOGbh . 



In this paper, we present a study using the Canada-France-Hawaii Telescope 
(CFHT) Legacy Survey data to search for cometary activity in dynamically 
asteroidal bodies. The first (smaller) data set is subjected to stringent tests 
designed to detect weak activity. These objects are examined by comparing 
their full-width-half-maximum (FWHM) measurements and brightness pro- 
files to those of stars of similar magnitude. If the profile of the asteroid is 
broader than those of the stars, it may indicate the presence of a coma. For 
such objects, the star profiles are subtracted from the profile of the asteroid, 
and an unpaired T-test is used to compare the residuals with the background 
to determine whether they are significantly different. A second (larger) data 
set is visually examined for evidence of stronger cometary activity. Using re- 
sults from both methods, upper limits are derived for the number of weakly 
and strongly active MBCs and are compared to those expected from collisional 
activation of asteroids. Finally, we present measurements of the known MBC 



2 



176P/LINEAR in the CFHT archive. 



2 Observations and Data Reduction 



All images were acquired with MegaCam on the 3.6-m CFHT in Hawaii. The 
images were taken as part of the Very Wid e (VW) segment o f the CFHT 
Legacy Survey (CFHTLS; I Jones et al.l fl2006h : iKavelaars et~aD fj2009h ). The 
asteroids forming the first data set were observed on 2004 December 15 — 16, 
2005 January 16 - 17, and 2006 May 1 - 2 and 25 - 26. All observations 
were taken in either the g' or r' filters, with exposure times of 90 s and 110 s, 
respectively. The limiting magnitude for a 90% probability of a 3a detection 
in the g' filter was 22.5, and 21.75 for the r' filter. The average seeing size was 
1" in both filters. 



The observations for the second data set were taken on various dates between 
August 2003 and January 2008. The exposure times for the g' and r' filters 
were 70 — 110 s and 110 — 180 s, respectively, yielding comparable limiting 
magnitudes to the observations above. 



The i mages were pre-processed using the Elixer pipeline (jMagnier and Cuillandre 



20041 ) . A fine astrometric correction was applied by the TERAPIX data- 



processing center, and the data was stored at the Canadian Ast ronomical Data 
Centre (CADC). Source Extractor (IBertin and Arnoutd . Il996l ) and additional 
software designed for detecting moving bodies were used to find asteroidal 
objects. A more detailed explanation of th e observations a nd da ta reduction 
process for the smaller data set are found in lWiegert et al.l (120071 ). which used 
that data to investigate the size distribution of kilometer-sized main-belt as- 
teroids (MBAs) function of color. 



For both data sets, each field was typically observed three times on the first 
night, approximately 45 minutes apart, and once the following night. This 
allowed for a reasonable determination of orbit parameters for main-belt ob- 
jects. The observations were taken at opposition and were not selected based 
on their possible main-belt content. Unlike other searches for MBCs, which 
focus on individual objects (Chamberlin et al. . 1996 ; Luu and Jewitt . 1990l ) 



or asteroid families (IHsieh and Jewittl . l2006bl ). this was a relatively unbiased 
survey of the main-belt. This allowed for a determination of the upper limit 
of MBCs expected for the whole main-belt, rather than for a sub-population. 
In this study, observations from the first night were utilized to search for 
cometary activity. When a second night of data was acquired, it was used for 
orbit refinement. 
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3 Data Analysis 



3. 1 Three-Level Analysis 



To determine whether an asteroid showed evidence of cometary activity, three 
levels of analysis were chosen, each refining the number of possible MBC can- 
didates. The investigation began with 1468 asteroids that were determined 
to be small bodies i n the main-belt by visual inspection of image triplets 



( jWiegert et all 120071 ) (each of these objects were also included in the visual 
investigation discussed below). The analysis was limited to objects with mag- 
nitudes less than 21.5 since the wings of the seeing profiles of fainter asteroids 
became too noisy for reliable examination in subsequent stages. Applying this 
limit left 952 objects to be examined. 

In the first level of analysis the FWHM measurements of the objects and 
stars were compared. Each asteroid image was rotated by the angle of the 
direction of motion, calculated from the RA and Dec coordinates, such that 
the direction of motion lied on the horizontal axis. The FWHM was measured 
perpendicular to the direction of motion using IRAF's IMEXAMINE tool. 
The FWHM along the direction of motion of the asteroid was not measured, 
since the images were trailed. 

All stars in the image were divided into magnitude bins (0.25-mag wide) and 
the median FWHM of each bin was calculated. The FWHMs of each asteroid 
and stars in the same magnitude bin were compared. If the FWHM of the 
asteroid was greater than 110% of the FWHM of the stars in two or more 
observations the asteroid was passed to the next analysis level. The arbitrary 
limit of 110% was chosen to provide a substantial number of asteroids with 
the highest FWHMs. Of the initial 952 asteroids, 415 passed this test. 

In the second level of analysis, the brightness profile of each asteroid was 
compared to three stars of similar magnitude (±0.25 mag) on the same image 
and CCD chip. This eliminated the effect of imaging variances across the CCD. 
The brightness profiles were created from a 20 x 20 matrix of pixel counts. 
The counts of the columns were averaged, the background was subtracted and 
the columns were plotted as a function of the line number. Each profile was 
subsequently scaled to a height of unity. 

The profiles of the asteroids were examined for interesting features (i.e., ex- 
tended wings) when compared visually to the stars. Although this process 
was subjective, the intention was to select objects having scaled counts in the 
wings more than 0.05 higher than the stars. An example of a typical extended 
profile that passed this test is shown in Fig. 1. There were 65 asteroids that 
showed a similar feature to Fig. 1, and subsequently passed this test. 



4 



[Figure 1] 



In the third level of analysis, the three star profiles were subtracted from the 
asteroid profile and an unpaired T-test was used to compare the residuals with 
the background noise. Using this method, 14 objects had residuals that were 
significantly different (i.e., T > 2 for a = 0.05) from the background in at 
least two observations. 

The images of those 14 objects were examined more carefully. All of the aster- 
oids were either overlapping or near a background object, internal reflections, 
or artifacts due to bright stars, thus making the object appear more extended. 
Therefore, no asteroids passed the final analysis level, or could be considered 
good MBC candidates. 



3.2 Visual Investigation 



The larger dat a set consisted of a total of 11615 objects which were visually 
inspected (as in lHsieh and Jewitt (j2006bh ). Some of these objects were outside 



the main-belt and were excluded from this analysis, leaving 11438 objects 
to be examined. A single known comet was found among the non-main-belt 
detections: active Centaur 166P/NEAT (a = 13.88 AU, q = 8.56 AU, e = 
0.383, i = 15.36°) in a distinctly active phase. 

An unknown object, observed on 2007 September 14, was discovered to show 
cometary activity. Its on-sky motion was consistent with being an MBC, but 
the orbit was poorly constrained. No known comets were found in the same 
region using MPChecker and Jet Propulsion Laboratory (JPL) Horizons. A 
report to the Minor Planet Center (MPC) did not produce any response linking 
it to another known body, either cometary or asteroidal; therefore, these are 
concluded to be the discovery observations of this object, pending further 
information from the MPC. It remains unclear whether this object is an MBC. 

Three observations of this object were taken over a 79-minute arc on one night 
and a fourth was taken approximately 48 hours later. The apparent magnitude 
varied from 21.2 to 21.6 in the r' filter. This yielded an absolute magnitude 
of 16.2 or a diameter, D, of 1.5 — 3.5 km for albedos ranging from 0.05 to 
0.25 (assuming the main-belt orbit calculated below). A faint, but distinct, 
coma/tail of approximately 10" in length was visible in all images (Fig. 2). 

[Figure 2] 

The images were retrieved from the CFHTLS archive at the CADC and 
searched for moving objects. Since this process was completed 12 months 
following the observations, telescopic follow-up at the time of discovery was 
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impractical. A search of the CFHT archive found no other images that were 
expected to contain the comet. 

The observations were taken near opposition (phase angle of 4°). The motion 
over 48 hours was consistent with an MBA and, assuming the usual Vaisala 
circular orbit, parameters of a ~ 3.9 AU and i ~ 13° were derived. However, 
the arc was too short to conclude this definitively, and the motion could also 
be fit with slightly smaller residuals (RMS 0.04" versus 0.14") by a more 
traditional cometary orbit with a ~ 6 AU, e ~ 0.6 and q ~ 2.5 AU. The 
nominal best fit cometary orbit would have passed perihelion in November 
2006, about a year before the object was observed in September 2007. 



Statistically, we expect this object is more likely to be a true comet rather tha n 
an MBC because of the relative rarity of the latter (IHsieh and Jewittl . l2006bl ). 
Its true nature will only be revealed once it is recovered. The observations, as 
reported to the MPC, appear in Table 1. 



[Table 1] 



The unknown comet and 166P were used to verify the three-level analysis tech- 
nique. Figures 3 and 4 show the brightness profiles of 166P and the unknown 
comet compared to three stars, respectively. Both objects passed the three 
levels of analysis. No other objects in this data set show cometary activity. 



[Figure 3] 
[Figure 4] 



4 Results 



4-1 Expected Number of Active MBCs 



Two of the three known MBCs (133P/Elst-Pizzaro and 176P/LINEAR) are 
members of the Themis family of asteroids in t he outer asteroid belt (3.05 < 



a < 3.22 AU, 0.12 < e < 0.19, 0.69° < i < 2.23°; IZappala et all fll990h ). while 



P/2005 Ul is just outside this family because of its slightly higher eccentricity 
(e = 0.25). The unknown comet from this survey resides in the outer belt 
(a ~ 3.9 AU), though its inclination of 13° is incompatible with the Themis 
family. 

Inside 3.2 AU, inter ior temperatures of astero ids are higher and may prevent 
the survival of ice ( IHsieh and Jewittl . l2006bl ). Therefore, when deriving the 
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expected numbers of MBCs, we initially assume that MBCs must be located 
in the outer belt, beyond 3.0 AU. 



The objects in this study range in size down to D ~ 1 km and D ~ 1.5 km for 
the large an d small data s ets, respectively. From the cumulative size distribu- 
tion given in lChengi (120041 ). the number of asteroids with D > 1 km is ~ 1 x 10 6 , 
while for D > 1.5 km is ~ 6 x 10 5 . Usi ng the de-biased num b er of asteroids 
as a function of magnitude derived by IJedicke and Metcalfe! (119981 ). ~ 30% 
of asteroids with absolute magnitude less than 16 have 3.0 < a < 5.0 AU. 
Extrapolating to smaller sizes, we estimate ~ 30% of all asteroids (~ 3 x 10 5 
with D > 1 km and ~ 1.8 x 10 5 with D > 1.5 km) are in the outer belt. 



In the first data set, 216 of 952 objects fall within 3.0 < a < 5.0 AU, which 
is consistent with the 30% figure stated above. Assuming MBCs only occur 
in the outer belt gives an upper limit of 1/216 (0.46%). Using the number of 
asteroids expected in this region (1.8 x 10 5 ) gives an upper limit of 830 ± 86 
currently weakly active MBCs, where the uncertainty is ±3cr (assuming a 
Poisson distribution). 



Using the 11438 objects from the second data set, the limit is constrained 
further. Although visual inspection may be insufficient to detect faint comae 
associated with weak activity, it is capable of determining stronger activity, 
as demonstrated by the detection of two moderately active comets. Of the 
11438 objects, 2667 have orbits with 3.0 < a < 5.0 AU, giving an upper limit 
of 0.037%, or 110 ± 31 str ongly active MBCs. This limit agrees with the the 
15 — 150 range derived by lHsieh and Jewittl (l2006bl ). 



If MBCs can occur anywhere in the main-belt, the active fractions for weakly 
and strongly active MBCs change to 0.11% and 0.009%, respectively. Using the 
total number of 6 x 10 5 for asteroids in the whole main-belt with D > 1.5 km, 
the upper limit for weakly active MBCs becomes 630 ± 77. For strongly active 
MBCs, the upper limit is 87 ± 28, assuming there are 1 x 10 6 asteroids with 
D > 1 km. 



4-2 Comparison with Collisional Lifetimes 



The upper limits derived for MBCs are compared with those expected from 
collisional activation of MBAs. The activation timescale for MBCs is unknown, 
but if collisional, is most likely comparable to the time between significant sub- 
catastrophic impacts . This r ate is approximately the destructional lifetime, 
which is derived by IChengl (120041 ) to be 10 8 years for asteroids with D ~ 
1 — 5 km (the size range of the majority of our objects). This is chosen as a 
first estimate of the collisional activation interval of MBCs. 
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Assuming there are 3 x 10 5 asteroids in the outer main-belt with D ~ 1 — 
5 km, and the act i ve lifet i mes of these objects after a collision is 10 3-4 years 



(IHsieh and Jewittl . l2006bl ; iDfaz and Gil-Hutton! 120081 ) . there should be 3 



30 active asteroids due to collisions. This is consistent with the upper limit 
derived for active MBCs. This derivation assumes MBCs remain consistently 
active for thousands of year s after activa tion, though observations indicate 
their activity is intermittent (jHsiehl (120071 ); see below). 



4.3 Mam-Belt Comet 17 6P /LINEAR 



The CFHT image archive was searched for observations of the three known 
MBCs. Only one was found, 176P/LINEAR, which was observed on 2005 June 
10 and 2007 January 15. For analysis, three 60-s observations from 2005 June 
10 were used, with seven and two minutes between observations, respectively. 
There were three 150-s observations from 2007 January 15, taken ~ 40 minutes 
apart. 

The profiles of 176P and stars of similar magnitude were compared, and the 
unpaired T-test was performed on both sets of observations. None of the im- 
ages passed the T-test, indicating 176P was not active during these times. 



Hsiehl (120071 ) determined 176P was active during November and December 
2005, just after perihelion. However, it was observed to be inactive shortly 
before and after this, in October 2005 and February 2006. With an orbital 
period of 5.71 years, 176P would be expected to be inactive during June 2005, 
when it was approaching perihelion, and January 2008, when it was about half 
way between aphelion and perihelion. Figure 5 shows the orbital position of 
the active and inactive phases of 176P. 

[Figure 5] 



5 Conclusions 



In this study, 952 asteroids are examined using three levels of analysis and 
11438 are inspected visually. In the larger set, one unknown object is discov- 
ered to display cometary activity. Since the orbit determination for this object 
has large uncertainties, it is possible it is an MBC. Statistically, however, it is 
more likely to be a true comet. The unknown comet, as well as 166P, are used 
to successfully test the three-level analysis technique. No other object shows 
evidence of cometary activity, so there are zero or one MBC candidates in the 
survey. 
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An upper limit is derived for weakly active MBCs to be 830 ± 86 for objects 
with D > 1.5 km and 3.0 < a < 5.0 AU. The larger sample of visually 
inspected objects in the outer main-belt with D > 1 km produces an upper 
limit of 110 ± 31 strongly active MBCs. Expanding to the whole main-belt, 
upper limits are 630 ± 77 and 87 ± 28 for weakly and stron gly active MBCs, 



respect ively. These agree with previous limits derived by iHsieh and Jewitt 



( I2006bl ) and are in good agreement with what is expected from collisional 
activation of asteroids. It is concluded that the activity of MBCs is consistent 
with activation caused by collisions with other small main-belt objects. 
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Observations of the Unknown Comet 



Obs. Date (UT) 


RA 




Dec 






mag 


2007 09 14.27098 


22 


51 


27.56 


+09 


59 


35.9 


21.3 


2007 09 14.30181 


22 


51 


26.38 


+09 


59 


29.6 


21.2 


2007 09 14.32586 


22 


51 


25.45 


+09 


59 


24.8 


21.2 


2007 09 16.29386 


22 


50 


12.88 


+09 


52 


33.4 


21.6 



Table 1 

Four observations of the unknown comet discovered in the CFHTLS data. Each 
image was taken in the r' filter and had exposure times of 180 s. 
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Line Number 

Fig. 1. The solid black line represents the brightness profile of an asteroid in the 
sample, while the three dotted lines represent the profiles of three comparison stars 
on the same chip. The asteroid shows slightly more counts in the wings, possibly- 
due to a coma. 
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Fig. 2. A combination of three images, taken 45 minutes apart, of the unknown 
comet observed 2007 September 14. The comet moves East to West in the sky. A 
tail of ~ 10" is seen to the right (West) of the object. See online supplementary 
material for an animation of this object. 
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Fig. 4. Brightness profile of the unknown comet (solid line) and three comparison 
stars (dotted lines). 
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Fig. 5. Orbital diagram of 176P showing its inactive and active phases. Orbits of 
Mercury, Venus, Earth, Mars, and Jupiter are shown in black, while the orbit of 
176P is in gray. Aphelion (A) and perihelion (P) positions are also marked. Solid 
and open squares show positions where 176P is active and inactive, respectively. 
The two positions marked with da tes are from the current study; all other positions 
are from lHsieh and Jewittl (|2006bh or iHsiehl (|2007h . 
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